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Abstract. The Tayler instability is a kink-type, current driven instability that plays 
an important role in plasma physics but might also be relevant in liquid metal 
applications with high electrical currents. In the framework of the Tayler-Spruit 
dynamo model of stellar magnetic field generation [T], the question of spontaneous 
helical (chiral) symmetry breaking during the saturation of the Tayler instability 
has received considerable interest [U El H. Focusing on fluids with low magnetic 
Prandtl numbers, for which the quasistatic approximation can be applied, we utilize 
an integro-differential equation approach in order to investigate the saturation 
mechanism of the Tayler instability. Both the exponential growth phase and the 
saturated phase are analyzed in terms of the action of the a and /3 effects of mean-field 
magnetohydrodynamics. In the exponential growth phase we always find a spontaneous 
chiral symmetry breaking which, however, disappears in the saturated phase. For 
higher degrees of supercriticality, we observe helicity oscillations in the saturated 
regime. For Lundquist numbers in the order of one we also obtain chiral symmetry 
breaking of the saturated magnetic held. 
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1. Introduction 

Electrical currents through an incompressible, viscous and resistive liquid conductor 
produce azimuthal magnetic helds which, beyond a critical held strength, become 
unstable to a non-axisymmetric, i.e. kink-type instability that we will call Tayler 
instability (TI) as a tribute to the seminal contributions of R.J. Tayler [El [7]. 

For a constant current density in an inhnitely long cylinder Rudiger et al had shown 
m that the governing parameter is the Hartmann number, Ha = B^(R)R{a/puy^‘^, 
which has to exceed a value in the order of 20 for the TI to set in {B^{R) is the azimuthal 
held at the outer radius R of the cylinder, a, p and v are the conductivity, density and 
viscosity of the huid, respectively). This critical value of Ha is actually consistent with 
previous results [101 HB [IB EE] concerning the ehects of viscosity and resistivity on 
the stability of various plasma z-pinches, if one leaves aside the ehects of complicated 
boundary conditions and non-homogeneous material parameters in the plasma case. 
Note that at these early days na it was far from obvious that the governing parameter 
for the onset of the instability is Ha, rather than the Lundquist number S = HaPm^/"^ 
(with Pm = vpoa denoting the magnetic Prandtl number, where po is the magnetic 
permeability constant). 

Whilst the focus of fusion related pinch experiments was prominently on the plasma 
destabilization when the ratio of axial to azimuthal magnetic held (the so-called safety 
parameter) falls below a certain critical value [H], a recent liquid metal experiment, 
with uniform conductivity and viscosity as well as well-dehned insulating boundary 
condition, has indeed conhrmed the Tl-threshold of Ha ~ 20 [T5|. From the application 
point of view, current-driven instabilities in liquid metals are presently considered a 
possible limitation for the integrity of large-scale liquid metal batteries. Such batteries 
are self-assembling stratihed systems made of a heavy liquid metal or metalloid (e.g., 
Bi, Sb) at the bottom, a suitable molten salt mixture as electrolyte in the middle, and 
a light alkaline or earth alkaline metal (e.g., Na, Mg) at the top. While small versions 
of liquid metal batteries have already been tested [lEl EZl EH EH], the occurrence of 
the TI could possibly present a serious problem for the stratihcation in larger batteries 
with prospected charging/discharging currents of some thousand amps. In [20] we had 
advised a simple trick to suppress the TI in liquid metal batteries by just returning 
the battery current through a bore in the centre. By the resulting change of the radial 
profile B^{r) it is possible to avoid the (ideal) condition d{rB‘^{r))/dr > 0 [7] for the 
onset of the TI. 

In a follow-up paper [5], a numerical code has been presented that is capable of 
treating Tl-problems at small values of Pm as they are typical for liquid metals. This 
was achieved by replacing the solution of the induction equation for the magnetic held 
by applying the so-called quasistatic approximation ED- This approximation allows to 
avoid the explicit time stepping of the magnetic held by computing the electrostatic 
potential by a Poisson equation, and deriving from this the electric current density. The 
induced magnetic held is then computed from the induced current density via Biot- 


The Tayler instability at low magnetic Prandtl numbers 


3 


Savart’s law. This way one arrives at an integro-differential eqnation approach, as it 
had already been used by Meir and Schmidt for different magnetohydrodynamic (MHD) 
problems Our numerical scheme utilizes the open source CFD library OpenFOAM® 
[23], supplemented by an MPI-parallelized implementation of Biot-Savart’s law. This 
code was then applied to a number of TI related problems, in particular for determining 
the scaling properties of the growth rate and the saturated velocity held, the dependence 
of the critical current on the geometric aspect ratio, as well as for validating various 
methods of preventing TI in liquid metal batteries m- Recently, our results were 
conhrmed by another code working completely in the framework of the differential 
equation approach, by analyzing the scaling properties of the solutions with Pm 125 ]. 
The authors also discussed carefully the limitations of the quasistatic approach for higher 
values of Pm. 

An interesting by-product of the battery-oriented simulations |5] was the 
observation of the transient occurrence, but ultimate disappearance, of helical structures 
during the evolution of the TI. On the hrst glance, the appearance of helical structures 
is surprising, since the underlying equations have no preference for left or right handed 
solutions. Yet, it is exactly this helical (or chiral) symmetry breaking that has gained 
considerable interest in various astrophysical problems. This applies in particular to the 
concept of the Tayler-Spruit dynamo [1] in which an azimuthal magnetic held is thought 
to become strong enough to drive the TI against the stable stratihcation in the radiation 
zone of a star. Combined with the usual differential rotation this effect might lead to 
a working dynamo. Despite of the attractiveness of the TI, in particular for explaining 
angular momentum transport in various types of stars [911261127I, the concept of the 
actual Tayler-Spruit dynamo is not without caveats. Zahn et al [2] have argued that 
the Tl-produced non-axisymmetric (m = 1) poloidal magnetic held alone would not be 
suited to close the dynamo loop (since the toroidal held wound up from it would have 
the same m = 1 dependence), but that some sufficiently large mean-held a ehect would 
be needed to produce the necessary axisymmetric poloidal held. 

It is exactly here where the question whether TI saturates with a hnite helicity, 
produced by a hnite a ehect, becomes highly relevant. Some recent papers have answered 
this question affirmatively: Gellert et al [3] have found spontaneous chiral symmetry 
breaking of the TI in simulations with Pm of 0.1, 1, and 10. Bonanno et al got a 
similar result for very large Pm = 10^. In addition to the numerical simulation, the 
latter authors developed a simple model of energy and helicity evolution resulting in an 
instructive phase portrait. The equations describing this behaviour can also be linked 
to a similar chiral symmetry breaking in biochemistry where it refers to the selection 
of one of two possible forms of bio-molecules (mainly sugars and amino acids) that are 
mirror images of each other [28] . 

With this background, the main motivation for the present paper is the discrepancy 
between the simulations of [3llH|2^ and the preliminary result of our low Pm simulations 
[5] showing that helicity starts to grow but ultimately decays to zero. Given the different 
Pm at which the respective simulations were done, it is worthwhile to understand in 
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detail the saturation mechanism of TI in dependence on Pm. 

Actually, helical states have a long history in plasma physics, tracing back to 
the early work of Lundquist on ’’Magneto-hydrostatic helds” Specializing general 
pressure-balanced helds to force-free helds that satisfy (V x B) X B = 0, he found, hrst, 
that helds with V x B = a{r)B fulhll this demand, and second, that a(r) = const 
must be requested for the held to remain force-free during its time-evolution. For 
cylindrical geometry, Lundquist found that the force-free condition, i.e. the demand 
that the current is parallel to the held, is guaranteed by Bessel function prohles 
= AjQ^ar), = AJi{ar) (interestingly, the very same prohles for the velocity held 
turned later out to provide the most efficient dynamo of the Ponomarenko or Riga type 

EB). 

Soon after Lundquist’s work, Chandrasekhar and Woltjer |32] interpreted this 
Bessel functions solution in terms of achieving ’’maximum magnetic energy for a given 
mean-squared current density” or, alternatively, as a ’’state of minimum dissipation for 
a given magnetic energy”. Since Bessel functions also maximize the magnetic helicity 
for given magnetic energy (and magnetic helicity is a better conserved quantity than 
energy) a surge of work was devoted to understand how solutions of this kind can be 
achieved dynamically. This goes mainly under the notion of Taylor relaxation [33] , and 
has found great interest in connection with the reversed held pinch. Quite a number 
of workers have tried to understand Taylor relaxation from diherent thermodynamic 
principles, such as minimum entropy production [3l| or minimum dissipation rates 

[asiiaHiETiEHiianiiiolii]. 

One of the hrst applications of a general thermodynamic principle to plasma 
relaxation goes back to a note of Max Steenbeck [42] (relying, in turn, on an earlier idea 
of Compton and Morse H). Steenbeck’s principle states that in real gas discharges 
at hxed current the heat power, and thus the voltage drop between the electrodes, 
is minimized (somewhat surprisingly, this minimum-dissipation principle corresponds 
perfectly with the maximum entropy production rate principle [H] if one considers the 
total system including the current-stabilizing external resistor Ha.) 

Interestingly, it was also Steenbeck who was later to create the theoretical 
framework that nowadays allows for a deeper dynamical understanding of those 
somewhat vague thermodynamic principles. Mean-held magnetohydrodynamics (MHD) 
was originally developed to explain self-excitation of cosmic magnetic helds [16]. Its 
main idea is that certain correlations of the small-scale parts of velocity and magnetic 
held contribute to the dynamics of the large scale magnetic held m- For helical 
turbulence the authors introduced the celebrated a ehect which drives an electromotive 
force parallel to a prevailing large scale magnetic held B. Similarly, turbulence leads to 
an increase of the resistivity by the /3 ehect, so that the mean electromotive force can 
be written in the form S = aB — fdV x B. 

Nowadays, mean-held concepts play not only a role in dynamo theory but also in 
the description of magnetically driven instabilities. Flow-driven helical dynamos and 
magnetically dominated helical “dynamos” are presently considered as two diherent 
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aspects of the very same mean-field MHD |3Sj. The detailed saturation mechanism of 
the TI, in particular its termination in a helical or non-helical state, is but one interesting 
application of mean-held MHD. 

In this paper, we are going to study the exponential growth and the hnal saturation 
of the TI in hnite cylindrical geometry for varying values of Ha and Pm. On the basis of 
an axisymmetric (m = 0) base state with an homogeneous axial current Jq that produces 
an azimuthal magnetic held Bq, we compute the m = 1 Tl-eigenmode comprising the 
velocity u and the induced magnetic held b from which we infer the mean electromotive 
force u X b (the overbar denotes the average over the azimuthal angle), and from this 
the mean-held coefficients a and /9. As a product of two m = 1 modes, u x b comprises 
certain m = 0 components that drive an azimuthal current (by virtue of the a ehect) 
and reduce the impressed axial current (by virtue of the (3 ehect). Although there is 
not a big scale separation between the m = 0 base state and the m = 1 perturbation, 
mean-held theory perfectly applies here. The electromotive force £ in direction of the 
mean held B (i.e. the a ehect), will be interpreted in terms of its relation to the small- 
scale current helicity, £ ■ B = —jb/a + eb, that had been derived and utilized by 
diherent authors [iniEniEIlESlEB]. 

In case of S' 1, the modihed currents and helds could be expected to resemble 
the typical Bessel function structure as typical for Taylor relaxation. In this sense, the 
mean axial held produced by the a ehect would follow from the principle of minimum 
dissipation [36] . 

However, this type of saturation mechanism, which relies on changing - by mean- 
held induction ehects - the electromagnetic base state in such a way that it becomes 
just marginally stable against TI, does not apply for S' 1. In this case the magnetic 
Reynolds number Rm of the Tl-produced how is much too small to induce any signihcant 
changes of the original applied magnetic held. The saturation must instead rely on a 
modihcation of the hydrodynamic base state, which we will discuss in detail. We will 
also evidence the occurrence of helicity oscillations, whose amplitudes and frequencies 
in dependence on Ha and Pm we will characterize. 

The paper closes with a discussion of the results, and with an outlook towards an 
application to stellar dynamo theory. 

2. The numerical scheme 

The usual numerical schemes for the simulation of TI, which solve the Navier-Stokes 
equation for the velocity and the induction equation for the magnetic held, are working 
typically only for values of Pm down to 10“^, although in a recent work by Herreman 
et al this limit has been challenged [25]. 

Here, we circumvent the usual Pm limitations of these codes by replacing 
the solution of the induction equation for the magnetic held by invoking the so- 
called quasistatic approximation [21]. We replace the explicit time stepping of the 
magnetic held by computing the electrostatic potential by a Poisson equation, and 
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deriving the electric current density. However, in contrast to many other inductionless 
approximations in which this procedure is sufficient, in our case we cannot avoid to 
compute the induced magnetic held, too. The reason for that is the presence of an 
externally applied electrical current in the huid. Computing the Lorentz force term it 
turns out that the product of the applied current with the induced held is of the same 
order as the product of the magnetic held (due to the applied current) with the induced 
current. Here, we compute the induced magnetic held from the induced current density 
by means of Biot-Savart’s law. This way we arrive at an integro-diherential equation 
approach, as it had already been used by Meir and Schmidt [22] . 

In detail, the numerical model as developed by Weber et al [5] works as follows: it 
solves the Navier-Stokes equations (NSE) for incompressible huids 

f 

u + (u ■ V) u =—Vp + uAu-\—- and V ■ w = 0, (1) 

P 

with u denoting the velocity, p the (modihed) pressure, fi, = JxB the electromagnetic 
Lorentz force density, J the total current density and B the total magnetic held. The 
NSE is solved using the PISO algorithm and applying no slip boundary conditions at 
the walls. 

Ohm’s law in moving conductors 

j = cr (—V(p + u X B) (2) 

allows to compute the induced current j by previously solving a Poisson equation for 
the perturbed electric potential ip = (j) — Joz/a: 

A(p = V ■ (n X B). (3) 


In the following, we will concentrate on cylindrical geometries with an axially applied 
current. Then, after subtracting the (constant) potential part Joz/a, with 2 ; as 
coordinate along the cylinder axis, we use the simple boundary condition <^9 = 0 on 
top and bottom and n ■ V(p = 0 on the mantle of the cylinder, with n as the surface 
normal vector. 

The induced magnetic held can then be calculated by Biot-Savart’s law 




dV‘ 


, 3 {r') x{r- r') 


\r — r’ 


( 4 ) 


Since this is a costly procedure, we modify here the method of [5] slightly. Actually, 
equation (4) is applied only at the boundary of the cylinder, while the magnetic held in 
the bulk is computed by solving the vectorial Poisson equation Ah = x {u x B) 

which results from the full time-dependent induction equation in the quasi-stationary 
approximation. 

Knowing b and j we compute the Lorentz force /l for the next iteration. A how 
chart of this numerical procedure is shown in hgure [H For more details about the 
numerical scheme, see section 2 and 3 of |5]. 
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^ end ^ 


Figure 1. Flow chart of the simulation model, slightly modified with respect to that 

of [5]. 


3. Results 

For the sake of simplicity, we consider a cylindrical electrically condncting flnid with a 
ratio of height L to diameter 2R of 1.25. A cnrrent is applied from top to bottom, jnst 
by setting the electrical potential to constant (bnt different) valnes at the two faces. 
Note that we refrain from taking into acconnt any cnrrents in the electrodes at top and 
bottom which has been shown to lead only to minor modihcations of the resnlts [5^ . 
The side walls of the cylinder are considered as electrically insnlating. No-slip bonndary 
conditions apply to the velocity at all bonndaries. 

In the following, we will focns on three different cases. The magnetic Prandtl 
nnmber for the hrst two rnns is Pm = 10“®. Differing in the Hartmann nnmber (60 and 
100), these two rnns will show a qnite different behavionr of the helicity. Whereas for 
Ha = 60 the initially growing helicity nltimately goes to zero, for Ha = 100 we observe 
helicity oscillation in the satnration regime. For the third case, with the mnch higher 
Pm = 10“^ and Ha = 100, we will hnd a hnite and non-zero valne of the hnal helicity. 
At the end of this section, we will snmmarize the different ways of satnration. 























The Tayler instability at low magnetic Prandtl numbers 
3.1. Saturation with zero helicity 


Here, we choose Ha = 60 and Pm = 10 ® which results in a Lundquist number S = 0.06 
that is dehnitely low enough for applying the quasistatic approximation. 






0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 




Figure 2. Time evolution of various quantities for Pm = 10“® and Ha = 60. (a) 

- Reynolds number, (b) - Normalized Rms value of the induced magnetic field, (c) 

- normalized kinetic helicity, (d) - Relation of electromotive force and small-scale 
current helicity, (e) - Normalized mean a effect, (f) - Normalized mean /3 effect, (g) - 
Normalized mean axial field, (h) - Normalized mean axial current. 

Figure |2] exhibits the time dependence of various quantities that characterize the 
instability. The indicated dimensionless time tn is t normalized by the viscous time scale 
/v. Figure 2a, to start with, shows the evolution of the averaged Reynolds number 
of the flow arising from the initial state at rest: {Re) = (here, (...) denotes 

an average over the total volume rather than only over the azimuthal angle, which will 
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always be indicated by an overbar). We chose a logarithmic scale in order to evidence 
the exponential growth of the TI. Approximately at tn = 0.27, the instability starts 
to saturate. We have added the respective energy contents of the various azimuthal 
wavenumbers m = 0,1, 2. Evidently, the m = 1 mode is the dominating one throughout 
the evolution. However, approximately at = 0.17, both the m = 0 and m = 2 modes 
start to increase with the double growth rate as the m = 1 mode. These even modes 
result from the non-linear term of the NSE. Saturation sets in when the m = 0 and 
m = 2 modes have acquired an amplitude comparably to that of the m = 1 mode 
which ultimately brings the growth rate of the TI to zero. The corresponding evolution 
of the averaged induced magnetic held is depicted in hgure 2b. Note that the m = 0 
component is here signihcantly weaker than the m = 2 component, quite in contrast to 
the rather parallel evolution for the kinetic energy. 

The kinetic helicity Hu = u ■ (V x u) is the next quantity to be discussed (see 
hgure 2c). Actually, we show here the helicity as normalized by the mean square of the 
velocity, i.e. {Hu)^ = {u ■ (V x u))R/{u^). After an initial increase, this normalized 
helicity stays nearly constant for a while (i.e., the non-normalized helicity {H^) grows 
with the same growth rate as the kinetic energy), until it decays to zero when saturation 
is reached. 

In Figure 2d we give an interpretation of the mean electromotive force 8 in 
direction of the large scale magnetic held in terms of the small-scale current helicity 
j ■ b, according to the relation S B = —jb/a + eb |l9l [50l EH EH ES]- We see 
that this relation is perfectly fulhlled, showing that a is essentially proportional to the 
current helicity, with a minor correction coming from an electric held term (note that 
all quantities are normalized here by BoJo/a). 

Figure 2e depicts separately the a ehect, dehned hj a = {u x b) ■ Bq/Bq. Since a 
has the dimension of a velocity, we give it here in form of a magnetic Reynolds number 
which includes again a complete spatial average: (a)^ = fioaR{{u x b) ■ Bq)/B^. We 
observe initially an (exponential) increase, though to a very small value in the order of 
and then a decay to zero. 

More monotonic than a, is the time evolution oi (3 = {u x b) ■ Jq/{^qJq) ehect 
(hgure 2f), which we normalize here by the magnetic dihusivity (/iocr)“^. This 
normalized, and spatially averaged (/?)„ = yiocr{{u x b) ■ Jq) / acquires values of 
about 6 X 10“^°, so its inhuence on the total resistivity can be considered as negligible. 

The induction ehects of the mean held coefficients a and (3 are illustrated in hgure 
2g and 2h. Figure 2g shows the mean axial magnetic held (6^) which is produced by the 
azimuthal current that is driven, in turn, by a. Normalized to Rq, we see again that this 
induction is negligibly small. Note also that (6^), in contrast to a (hgure 2e), does not 
completely vanish in the saturation regime. We attribute this to numerical inaccuracies 
which seem to “stray” some energy from the much stronger m = 1 mode into the m = 0 
and m = 2 modes, an ehect that is already visible in the non-physical parallelism of all 
modes in the beginning of the exponential growth regime (hgure 2a). In contrast to this 
slight discrepancy, the behaviour of {jz)/Jo (hgure 2h) is nearly identical to that of /3 
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(figure 2f). 



b. 













Figure 3. Snapshots of different quantities for three different instants. The contours 
correspond to 7%, 3% and 50%, 20%, and 10% of the extremal values in the first, 
second, third and fourth row, respectively. 


For three selected instants in time, figure EliHustrates the spatial structure of various 
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features of the TI (the normalized values of u^-, {u x b) ■ Bq, {u x b) ■ Jq). The 
left column depicts these quantities amidst the exponential growth phase in which we 
observe a clear helical structure. The middle and right columns show then the respective 
structures shortly before and during the saturated regime. Evidently, the helicity has 
completely disappeared here. In terms of the mode structure, we notice that the left 
handed spiral (m = 1, say) and the right handed spiral (m = —1, say) have grown to 
the same strength. 

Having seen that, due to the low values of Pm, neither a nor /3 is able to induce 
any relevant change of the electromagnetic base state that would lead to saturation of 
the TI, we have to look for an alternative saturation mechanism. Evidently, this can 
only be related to a change of the hydrodynamic state, i.e. the flow held. Let us return 
to hgure [2^: after a long exponential growth period of the (more or less) pure m = 1 
mode, at tn = 0.17 the m = 0 and m = 2 modes start to grow due to the action of the 
nonlinear terms in the NSE. As shown in hgure IU the m = 0 part of the velocity in 
the saturation comprises two poloidal vortices pointing outward in the equatorial plane 
(in the dynamo community this how topology would be denoted by s2+ |55]). This 
axisymmetric state, together with the m = 2 component, changes the hydrodynamic 
base state for the TI so that it becomes just marginally stable. 


Sto Tn = 0 m = 1 Tn= 2 



Figure 4. The velocity field in the saturated state, including the three lowest 
azimuthal modes. 

With the simultaneous appearance of an m = 0 and m = 2 component, it is no 
surprise that the saturation is also connected with a restoration of the chiral symmetry. 
According to the sum rule for the nonlinear interaction the m = 2 will produce, from 
any dominant m = 1 mode, a corresponding m = — 1 mode, so that chiral symmetry is 
ultimately restored. 

3.2. Saturation with helicity oscillations 

We now increase the Hartmann number from Ha = 60 to Ha = 100. As before, 
hgure [5] illustrates the time evolution of various quantities. While the behaviour of 
the Reynolds number and the (3 effect are not signihcantly different from the previous 
case with Ha = 60, the helicity and a look quite differently. Evidently, the TI does 
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not anymore saturate with zero helicity, but instead produces a helicity oscillation in 
the hnal state. Figure 6 shows again the spatial structure of various quantities in the 
exponential growth phase and at two different instants in the saturated state. In the 
plots for {u X b) ■ Bq, the helicity oscillation appears as a slightly changing asymmetry 
between positive (orange) and negative (blue) “blobs”. 



-5x10-3 





Figure 5. Same as figure [21 but for Pm = 10 ® and Ha = 100. 

For the range between Ha = 40 and Ha = 140 we characterize this helicity 
oscillation by its amplitude and frequency (figure [7]). 

3.3. Saturation with finite helicity 

We leave now the realm in which the magnetic induction is irrelevant for saturation and 
move towards the parameter region which had already been explored by other codes. 
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Figure 6. Same as figure O but for Pm = 10 ® and Ha = 100. 


Actually, we increase the magnetic Prandtl number to Pm = 10“^, and consider the 
case Ha = 100. This leads to a Lundquist number of S' = 3.16 for which we come close 
to the edge of applicability of the inductionless approximation Eg. 
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Figure 7 . Characteristics of the helicity osciallations in dependence on Ha, for 
Pm = 10“®. (a) Frequency, also compared with the growth rate in the exponential 
growth phase, (b) Amplitude. 


Again, figures |8] and [9] show the time evolution, and some snapshots, of various 
quantities. The main difference to the former runs at low Pm is that now the helicity and 
a acquire non-zero values in the saturated state. This is also shown in the snapshots of 
£gure[9l Evidently, we are now in a regime where the induced magnetic helds contribute 
signihcantly to the saturation mechanism. This involves also that the m = 0 component 
of the magnetic held (I8 ]d) becomes now comparable to the m = 2 part, quite in contrast 
to the former cases with Pm = 10“®. 

This is illustrated in figure [10] which shows the dependence of the induced mean 
current {jz)/Jo and the induced mean axial held {b^) /Bq on Pm (at hxed Ha = 100). 
According to the criterion of Kruskal-Shafranov mi, we know that tends to inhibit 
the TI, so that a hnite value of a is indeed likely to appear in the saturation regime. 

3.4- Between chiral symmetry breaking and helicity oscillations 

In the following, we will summarize the three diherent saturation mechanisms. We will 
be guided by the simple and instructive model of chiral symmetry breaking model that 
had been worked out by Ron a n n o et al [1]. 

The authors started from left and right handed helical modes for the velocity and 
the magnetic held, fulhlling the Beltrami relation 

V X T = XL and V x i? = —XR (5) 

which can be realized by appropriate linear combinations of Chandrasekhar-Kendall 
functions cos(m0) cos{nz7i/H). 

Invoking some symmetry arguments, the authors “guessed” the simplest form of a 
Lagrangian which then leads to the following evolution equations for the energy of the 
left and right handed helical modes 

= 2'yEL — 4:^E‘l — A^^ElEr 

dEji 


dt 


‘2'yER — 4/iE^ — A^^^ElEr . 


( 6 ) 

( 7 ) 
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Figure 8. Same as figure [U but for Pm = 10 ^ and Ha = 100. 


The last terms on the r.h.s. of equations (6,7) describe the so-called mutual 
antagonism between the two chiralities that has been used extensively in the theory 
of homochirality of bio-molecules [56l |57l 1^ . 

From here, one arrives at the following evolution equations for the total energy 
E = Er -|- El and the helicity El = Er — Er- 

dE 

— = 2'yE - 2(/i -b p^)E‘^ - 2(/r - 

dt 

^ = 2-fH - AfiEH . 

In hgure [11] we illustrate the phase portrait of this equation system showing a 
clear chiral symmetry breaking both in the exponential growth phase as well as in the 


( 8 ) 

(9) 
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Figure 9. Same as figureO but for Pm = 10 ^ and Ha = 100, and different instants, 
saturated phase. 

We return now to our three cases, for which we plot the time evolution (Figure [121 
left columns) and the phase portrait of the kinetic helicity (middle column), as well as 
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Figure 10. Dependence of (&z)/-Bo and {jz)/Jo on the magnetic Prandtl number 
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Figure 11. Phase portrait of the coupled equation system (8,9), for the parameters 
7 = 2.71, y = 3.0, /r* = 5.7. In this typical situation. Si is a repeller, S 2 and S 3 are 
attractive points, while 54 is a saddle point. 


the phase portrait of the current helicity (right column). To make contact with hgure 
11, we have now chosen a different normalization so that both helicities start at zero. 

We start, in the hrst row of hgurewith Pm = 10“® and Ha = 60. Evidently the 
exponential growth phase looks very similar as in hgure 11, but ultimately the system 
runs into a state with zero helicity. 

The second case with Pm = 10“®, Ha = 100 is similar but terminates with a 
helicity oscillation around zero. It is remarkable that this helicity oscillation proceeds 
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Figure 12. Time evolution (left), and phase portraits for the kinetic helicity (middle), 
and phase portrait of the current helicity (right), for four different cases. Note that 
the viscous time scale is indeed the relevant one both for the growth process and for 
the helicity oscillations. 


without any significant oscillation of the energy. 

The third case, Pm = 10“^, Ha = 100 has indeed some resemblance with the above 
model of Bonanno et al and terminates with a finite, non-zero helicity. 

In the fourth row, we add here also the plot for Pm = 10“^ and Ha = 100 
which amounts to S' = 10 which is beyond the simple applicability of the quasistatic 
approximation and should be, therefore, considered with caution. Evidently, the 
magnetic helicity shows now a phase portrait that is similar to figure [TTl despite the 
fact that we now observe an ’’overshoot” to the helicity with the other sign, and also a 
remaining oscillation in the saturated state. While we do not claim that this is a reliable 
result, the two last rows of figure HU at least suggest that we are now approaching the 
usual saturation scheme as already discussed by |1]. 
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4. Conclusions 

In this paper, we have utilized an integro-differential equation solver for addressing 
the problem of chiral symmetry breaking in the exponential growth phase and in the 
saturation phase of the Tayler instability. The advantage of this code is its easy 
applicability for small magnetic Prandtl numbers, while its suitability for problems at 
S' > 1 is at least questionable (at least it has to be checked case by case whether a final 
steady state with finite and non-oscillatory, i.e. static, helicity could still be treated 
with our quasistatic scheme). Our simulation have allowed to identify three different 
saturation regimes. 

To start with the last regime, for a comparably large value of S' = 3.16 we have 
confirmed a similar type of chiral symmetry breaking as it was previously evidenced by 
Gellert et al [3] and Bonanno et al [1]. Depending on the random initial conditions, 
the TI grows with one of the two possible helicities which does not disappear in the 
saturated regime. The helicity is intrinsically connected with a non-zero a effect that 
generates a current parallel to the applied azimuthal magnetic field. At the same time 
the mean-field e.m.f. contains also a significant (3 effect that changes the axial current. 
Both effects together work against the TI. In the ultimate case of high S' (which is, 
probably, not accessible by our code) one could expect a sort of Taylor-relaxation into a 
helicity maximizing state [33]. Whether for those large values of S' one reaches a regime 
of helicity oscillation around a finite value (as suggested by the fourth row of figure [T2D 
is still to be validated by complementary codes. 

The saturation mechanism, which relies on changing (by mean-field induction 
effects) the electromagnetic base state in such a way that it becomes just marginally 
stable against TI, does not apply for S' -C 1. This can already be seen from the general 
scaling Re oc Ha^, or Rm oc S'^ which applies both to S' < 1 and S' > 1 [5|. For S' <C 1 
the final Rm becomes much too small in order to induce any significant changes of the 
original applied magnetic field. 

In this parameter region, the saturation relies instead on the non-linear appearance 
of an m = 0 and an m = 2 velocity component, which now changes the hydrodynamic 
base state of the TI in a way that the growth rate of the TI vanishes. 

Perhaps the most interesting result of our study is the observation of a robust 
and systematic helicity oscillation whose amplitude and frequency dependence on Ha 
has been worked out. Interestingly this helicity oscillation is not connected with any 
significant energy oscillations. 

Based on the latter observation, we would like to conclude with some, admittedly, 
very speculative ideas. There is a long tradition in trying to link the various frequencies 
of the solar dynamo action to corresponding periodicities of planetary motion. Tracing 
back to a paper by Jose [5H], who had related the 11.86 years Jupiter orbit with the 
22.08 year solar cycle, some refinements of this connection in terms of a combined torque 
and gravity action of Earth, Venus and Jupiter have been discussed recently [59]. Other 
papers have tried to link periodicities of the Jupiter-Saturn orbits to various longer- 
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time cycles of the solar dynamo [501 El], with possible connections to the climate on 
Earth. However, in all those cases, it was noticed that the planetary forces are mnch 
to weak to compete with the typical acceleration forces in the convection zone. The 
only viable way of influencing the solar dynamo was speculated to rely on the action 
of gravity on the shape, or local rotation rate, of the tachocline. Yet, this would imply 
that the solar dynamo works indeed as some sort of Tayler-Spruit dynamo |1], in which 
the transformation from poloidal to toroidal held is traditionally realised by differential 
rotation, the reverse mechanism, however, by some a effect due to the TI. 

It is exactly here where helicity oscillations, and their possible synchronization with 
planetary forces and torques, might come into play. In particular, since the oscillations 
of a are not connected to any signihcant changes of the energetic content, very minor 
changes of the state of the tachocline might just open the “a-bottleneck” for the Tayler- 
Spruit dynamo (which is, in any case, still a sort of a — H like dynamo). Even if 
the a oscillations are around some non-zero mean values prevailing in the two solar 
hemispheres, it still might give rise to dynamo oscillations. 

Note that a parametric resonance and synchronization of m = 1 dynamo 
eigenmodes with m = 2 velocity perturbations has been observed both for galactic 
dynamos (swing excitation, [62|), and for a VKS-like dynamo |63|. Whether a similar 
effect may actually be at work for synchronizing the solar dynamo with periodic 
planetary forces via their action on the tachoclinic state, will remain a topic for future 
investigations. 
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